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Abstract 

The stability of matter composed of electrons and static nuclei is 
investigated for a relativistic dynamics for the electrons given by a 
suitably projected Dirac operator and with Coulomb interactions. In 
addition there is an arbitrary classical magnetic field of finite energy. 
Despite the previously known facts that ordinary nonrelativistic mat- 
ter with magnetic fields, or relativistic matter without magnetic fields 
is already unstable when a, the fine structure constant, is too large it is 
noteworthy that the combination of the two is still stable provided the 
projection onto the positive energy states of the Dirac operator, which 
defines the electron, is chosen properly. A good choice is to include 
the magnetic field in the definition. A bad choice, which always leads 
to instability, is the usual one in which the positive energy states are 
defined by the free Dirac operator. Both assertions are proved here. 

1 Introduction 

The stability of matter concerns the many-electron and many-nucleus quan- 
tum mechanical problem and the question whether the ground state energy 
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is finite (stability of the first kind). If so, is it bounded below by a constant 
(which is independent of the position of the nuclei) times the number of par- 
ticles (stability of the second kind)? The linear lower bound is important for 
thermodynamics, which will not exist in the usual way without it. 

The first positive resolution of this problem for the nonrelativistic Schro- 
dinger equation was given by Dyson and Lenard |7], || and approached dif- 



ferently by Federbush [10|. The constant, i.e., the energy per particle, was 



considerably improved by Lieb and Thirring in J21], |22|. Following that, the 
stability of a relativistic version of the Schrodinger equation (in which p 2 is 
replaced by a/p 2 + m2 ) was proved by Conlon |J and later improved by Lieb 



and Yau [23| who showed that matter is stable in this model if and only if 
the fine structure constant a is small enough and if Zat <2/tt. (See [[£1 for 
a historical account up to 1995.) A recent result of Lieb, Loss, and Siedentop 



that we shall use is in flTjJ and is discussed in Section |[ 



In these works the nuclei are fixed in space because they are very massive 
and because we know that the nuclear motion is largely irrelevant for under- 
standing matter. In other words, if nuclear motion were the only thing that 
prevented the instability of matter then the world would look very different 
from what it does. We continue this practice here. 
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There is, however, a more important quantity that requires some atten- 
tion, namely magnetic fields. It was noted that the action of such fields on 
the translational degrees of freedom of the electrons p — > p + eA, can lower 
the energy only by an inconsequential amount. This is a kind of diamag- 
netic inequality. On the other hand, spin-magnetic field interaction (in which 
(p + eA) 2 is replaced by the Pauli operator [<x-(p + eA)] 2 = (p + eA) 2 + ecr-B 
can cause instability. The energy is then unbounded below if arbitrarily large 
fields are allowed, but this is so only because the energy of the magnetic field 
has not been taken into account. Does the field energy, (87r) _1 J B 2 , insure 



stability? This question was raised for the nonrelativistic case in and 
finally settled in a satisfactory manner in |2tJ (see also Bugliaro et al. || and 
Fefferman |TT|]). The upshot of this investigation is that stability (of both 
first and second kinds) requires a bound on both a and on Za 2 . 

Other related results are the stability of non-relativistic matter with a 
second quantized, ultra-violet cut-off photon field (Frohlich et al. |12| ). 

Both the passage to relativistic kinematics (which, in quantum mechanics, 
means that both the kinetic energy and the Coulomb potential scale with 
length in the same way, namely like an inverse length) and the introduction 
of the nonrelativistic Pauli operator require a bound on a and on Z for 
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stability. The combination of the two might be expected to lead to disaster. 
We find, however, that it does not necessarily do so! 

Our main result is that matter is indeed stable with a suitably defined 
relativistic kinematics. This is shown in Section || 

The proper way to introduce relativistic kinematics for spin- 1/2 particles 
is via the Dirac operator, but this is unbounded below. A resolution of this 
problem, due to Dirac, is to permit the electrons to live only in the positive 
energy subspace of the Dirac operator. This idea was further pursued by 
Brown and Ravenhall [[| (see also Bethe and Salpeter in their Handbuch 
article J!]) to give a quantitative description of real atoms. 

There are, however, other Dirac operators (which include electromagnetic 
potentials) whose positive subspace can be used to define the space in which 
the electrons can live. (To avoid confusion, let us note that the Hamiltonian is 
formally always the same and includes whatever fields happen to be present. 
The only point to be resolved is what part of the one-particle Hilbert space 



is allowed for electrons.) The review articles of Sucher [gj, [25|, |26| can be 



consulted here. These choices have also been used in quantum chemistry and 



other practical calculations, see, e.g., [14, 15 



All of these choices have in common that there is no creation of electron- 
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positron pairs explaining the name "no-pair Hamiltonian" for the resulting 
energy operator. (Note that we could also treat positrons or a combination 
of electrons and positrons interacting by Coulomb forces in a similar way.) 

There are three obvious choices to consider. One is the free Dirac oper- 
ator. This always leads to instability of the first kind when a magnetic field 
is added unless the particle number is held to some small value (see Section 
f|). Note also that this choice leads to a non-gauge invariant model: multi- 
plication of a state with the factor exp(i0(r)) for a non-constant gauge is not 
allowed, since it leads out of the positive spectral subspace. 

Remarkably, the Dirac operator that includes the magnetic field always 
gives stability, if Z and a are not too large, as in the two cases (relativistic 
without magnetic field and nonrelativistic with magnetic field) mentioned 
above (see Section ||). This model is gauge invariant. 

The third choice which, indeed, is sometimes used, is to include both the 
one-body attractive electric potential of the nuclei and the magnetic field 
in the definition of the Dirac operator that defines the positive subspace. 
If this is done then the question of stability is immediately solved because 
the remaining terms in the Hamiltonian are positive, and hence the total 
energy is ipso facto positive. This choice, which is important but trivial in 
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the context of this present inquiry, will not be mentioned further. 

Oddly, the instability proof given in Subsection |4.2| is much more compli- 
cated than the stability proof (Section ||). This reverses the usual situation. 



A preliminary version of this paper appeared in | 27fl ; the present version 



is to be regarded as the original one (as stated in J27j) and contains several 



significant corrections to the preliminary text in p7| . In particular, the proof 
of Theorem 2 and the first half of the proof of Theorem 1 have been corrected 
and simplified. 



2 Basic Definitions 

The energy of N relativistic electrons in the field of K nuclei with atomic 

numbers Zi,...,Zk e R + located at Ri,...,Rx € M 3 which are pairwise 

different in a magnetic field B = V x A in the state \I> is — following the 
ideas of Brown and Ravenhall || — 

N 

£ A [% ] := (¥, (J2 Du(A) + aV c )V) + — / B{vfdv. (1) 
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Here D V {A) := a ■ (—iV u + eA(r„)) +m/9 is the Dirac operator with vector 
potential A. Furthermore, 



N K N K 

v c -.= -j2Yl + + J2 

V=\ K = l ~ ^1 M- = l l 1 "^ _ ^ «,A=1 1^ _ R - 



(2) 



is the Coulomb interaction between the particles, and B(r) := |V x A(r)| 
is the modulus of the magnetic field. Planck's constant divided by 2n and 
the velocity of light, are taken to be one in suitable units. The fine structure 
constant a equals e 2 , where — e is the electron charge. Experimentally, a 
is about 1/137.037. The mass of the electron is denoted by m. The 4x4 
matrices ct and (5 are the four Dirac matrices in standard representation, 
namely 
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Finally, the state \& should have finite kinetic energy, i.e., it should be 
in the Sobolev space HV 2 [(R 3 x {1,2,3,4})*], and should also be in the 
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electronic Hilbert space of antisymmetric spinors 

N 

Hn,a '■= f \ , (3) 

where 7i + is the positive spectral subspace of the Dirac operator D(A) and 
where A is some vector potential to be chosen later. The vector potential A 
serves to define the positive subspace. Two choices will be considered here. 
One is A = 0, in which case we are talking about the free Dirac operator. 
This choice, or model, goes back to Brown and Ravenhall ||. As we shall 
see in Section |], the resulting energy functional — apart from being not gauge 
invariant — is not bounded from below. A natural modification of the model, 
namely to take A := A is not only gauge invariant, but will also turn out to 
be stable of the second kind (see Section D). 

The quantity of interest is the lowest possible energy 

En,k '■= inf £4 

where the infimum is taken over all allowed states all allowed vector 
potentials A, and over all pairwise different nuclear positions Ri, ...,R#. 

In the case of a single nucleus without a magnetic field, the energy form 
Sq was shown in || to be bounded from below, if and only if ctZ < ctZc ■— 
2/(7r/2 + 2/7r) > 2/7T , which corresponds to Z w 124. We will not be able 
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to reach this value in the general case of many nuclei and when the electron 
state space is not determined by the free Dirac operator. The reason is 
that special techniques were used in [f| to handle the one-nucleus case; these 
techniques took advantage of the weakening of the Coulomb singularities 
caused by the fact that states in 7i + cannot be localized in space arbitrarily 
sharply. Unfortunately, we do not know how to implement this observation 
with magnetic fields and many nuclei. 



3 Stability with the Modified Projector 



Our proof of the stability of matter when the vector potential A is included 
in the definition of the positive energy electron states will depend essentially 
on three inequalities: 

BKS inequality: For any self-adjoint operator X, the negative (positive) 
part, X T is defined to be (\X\ =p X)/2. Given two non-negative self- 
adjoint operators C and D such that (C 2 — D 2 )]/ 2 is trace class, we 
have the trace inequality 



This is a special case of a more general inequality of Birman, Koplienko, 




(4) 
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and Solomyak [TJ]; in particular, the number 2 in (|j) can be replaced by 
any p > 1. A proof for the special case of the inequality needed here is 
given in Appendix [A]. 

N 

Stability of relativistic matter: On /\ (if 1 / 2 (M 3 ) <g> O 3 ), the fermionic 
Hilbert space, we have 

N 

J2\~^u-A\+ aV c > 0, (5) 

(where | • • • | means \J (• • -) 2 ) for all vector fields A : M 3 — > M 3 with, 
e.g., square integrable gradients, if 

1/5 > 1/5 C := (n/2)Z + 2.2159g 1/3 Z 2/3 + 1.0307g 1/3 (6) 

and Z\, Z K < Z. 

We wish to use this inequality for 4-component spinors, i.e., q = 4. 
However, we are interested in the subspace TCn,a i n which the particles 
are restricted to the positive energy subspace of the Dirac operator, 
D(A). Although q — 4, the 'effective' q is really 2, and the analysis 
in Appendix is our justification for this. The only thing that really 
counts in deriving (||) is the bound on the reduced one-body density 
matrix 7 mentioned in Appendix [FJ. 
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The stability of the relativistic Hamiltonian (||) was first shown by 
Conlon H for A = 0. The best currently available constants with 



^4 = are in |23fl while (g), which is taken from |19[], is the best known 



result for general A. 

Semi-classical bound: Given a positive constant //, a real vector field A 
with, e..g., square integrable gradients, and a real- valued function <p G 
L 2 (1R 3 ) the inequality 

tr[(-i^-Ay-^<^L f <fr (7) 

holds, which is a special case of the Lieb-Thirring inequality (see | 22| , 
T7J). It is known that Li/2,3 < 0.06003. The left side of (0) is simply 
J2j l-^il 1 ^ 2 ) where the Xj are the negative eigenvalues of the operator 



As an illustration of the usefulness of the trace estimate @ , let us combine 
it with the Lieb-Thirring inequality (0) (or any other Lieb-Thirring inequality 
for that matter) to derive some previously known inequalities. The constants 
obtainable in this way are comparable with the best ones known so far: 

Daubechies inequality: We begin with a "relativistic" inequality that was 
first proven by Daubechies 0. By replacing by ip 2 in ([?]), we get using 
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tr (| - iV - A\ - cp)_ < Li/2,3 / (8) 

The constant 0.06003 obtained here should be compared with the num- 
ber 0.0258 in §. 

Non-relativistic magnetic stability: A non-relativistic analogue of our 
main problem is to bound the form 

N 

£ := (*, (V P„(A) + aV c )V) + — / B(r) 2 dr (9) 

which was treated in |20fl . Here P U (A) := [cr ■ (— iV u + eA(r i/ ))] 2 is the 
Pauli operator with vector potential A. 

First, we note that x 2 > +A|a;| — A 2 /4 holds. A constant in the energy 
form, however, is irrelevant for checking on stability of the second kind. 
Using (|) it is then enough to show the positivity of 

-tr(AP(A) 1/2 -/t| -iV + eA|)_ + — / B(r) 2 dv (10) 

where we have set k := a/a c . The trace in this and the next expression 
are over L 2 (R 3 ) ® C 2 . Using the BKS inequality gives the lower bound 

-tr[(X 2 -K 2 )\-tV + eA\ 2 )-e\ 2 B(r)} 1 / 2 + — [ B{v) 2 dv. 

87T J R 3 
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Applying the Lieb-Thirring inequality (|7]) yields the following sufficient 
condition for stability (recall that a = e 2 and that there are two spin 
states) 

Zljl o ; — \ . 

3. 3 (A2-k2)3/2 " 8 7l 

Optimizing in A gives 

16a 2 1 

which gives for the physical values a ~ 1/137.037 and q = 2 a range of 
stability up to Z < 1096, which is to be compared with Z < 1050 in 
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We turn now to our main result. 

Theorem 1 Let Zi,...,Zx < Z < 2/(ira) and let a < a c where a c is the 
unique solution of the equation 



(167rLi i3 « c ) 2 / 3 = l-a 2 /d c 2 , 



with a c := [(n/2)Z + 2.2159 • 2 1 / 3 Z 2 / 3 + 1.0307 • 2 1 / 3 ]" 1 as in g). Then £ A 
is non-negative. 

Numerically, this gives 

Z < 56 
14 



when evaluated with the experimental value a m 1/137.037 for the fine struc- 
ture constant. Alternatively, considering hydrogen only, i.e., Z — 1, we ob- 
tain the upper bound 

a < 1/8.139 

for the fine structure constant. It is a challenge to improve this result so 
that it covers all physical nuclear charges and the physical value of the fine 
structure constant, as was done for K — N — 1 and A = in ||. 

Proof. The first step in our proof is to utilize (Q) to replace V c by the 
one-body operator (— l/a c ) $^„=i I — ^Vj, + eA|, where a c is given by (^) 
with q = 2, as we explained just after (|G|). (The idea of using the relativistic 
stability result (^) to bound the Coulomb potential by a one-body operator 
first appears in |2(|.) Our energy £a is now bounded below by 

N 

) := (viz, V(A,(A) - k| - iV„ + eA(r„)|)*) + — / 5(r) 2 rfr, (11) 
where k := a/a c . 

The first term on the right side of (|TT| ) is bounded below by the sum of 
the negative eigenvalues, — trh_, of the one-body operator 

h ■= A+(d(A) - k\ - zV + eA(r)|)A+, 

where A + is the projector onto the positive spectral subspace of D(A). 
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Let us define 

S := \D(A)\ - k\ -iV + eA(r)|, 

whence h = A+ SA + , because A + L>(A)A+ = A + |D(A)|A + . We note that 
for any two self-adjoint operators X and p with X > and < p < 1, 
trX > trpX. With p taken to be the projector onto the negative spectral 
subspace of h we then have that 



trh- = —trph = —trpA + SA + 

= trpA + S-A+-trpA + S + A + (12) 
< trA +1 S_A+. 

We introduce the 4x4 unitary 




and note that IJ- 1 D{A)U = -D(A). Therefore, U^A+U = A_. 

It follows from the spectral theorem that for any self-adjoint X, unitary 
U, and function F 

FiU^XU) = U- l F(X)U. 
With F{t) = \t\, we then have that IJ- 1 \D(A)\U = \D(A)\, and hence 

u- l su = S. 
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Therefore, since U^A+U = A_, and with F(t) = \ - t) = t_, we 
have that U~ 1 S-U = S- and 

tr A + S_A + = tr A_S_A_. 

Hence, using (f[|), 

tr/i_ < itr(A + S_A + ) + - tr A_S_A_ = -tiS— 

_ & Zi 

(Note: much of the preceding discussion was needed only to get the factor 
1/2 here. This factor improves our final constants for stability.) 
Next, we use the BKS inequality (§) to bound trSL as follows: 



tr h- < - tr £>_ < - tr 
~ 2 ~ 2 



D(A) 2 - k 2 \ - zV + eA(r)| 2 ^ \ (13) 



Y 

However, D(A) 2 = | | with Y = -P(A) + m 2 , and where 

o y r 

P(A) = [<r • (-iV + eA)] 2 = | - iV + eA(r)| 2 + e<7 ■ B(r) 



is the Pauli operator. 

Since X i— > trXi^ 2 is operator monotone decreasing, we see that our 
lower bound for the energy is monotone increasing in m, and thus it suffices 
to prove the positivity of TCn,a in the massless case. The key observation is 
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that our lower bound involves only tr SL in the entire one-body space, not 
the positive energy subspace. The energy would not be obviously monotone 
in m if we had to restrict functions to the positive subspace, since changing 
m would also entail changing the space. This problem does not arise in the 
absence of the positive subspace constraint. 

Because of the 'diagonal' structure of the operator S, we can drop the 
factor 1/2 by replacing the trace on L 2 (IR 3 ) <g>C 4 by the trace on L 2 (M 3 ) <g> C 2 . 
This yields 



We regard the operator in the second line as acting on functions (of one 
component only) instead of spinors, which accounts for the factor two (and 
not one).. Finally we apply the Lieb-Thirring inequality ([?]) to the right hand 
side yielding (recall that e 2 = a) 





B{rfdY. 



(14) 



£a[*, *] > [-2^1,3 a(l - a 2 /a 2 r 3/2 + l/(8rr)] / B 2 {r) 2 dr. 



Jr 3 



Thus we need 



(167rLi 3 a) 2/3 < 1 - a 2 /a 2 . 



(15) 
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Since the right hand side of this inequality is monotone decreasing in a for 
positive a, while the left hand side is monotone increasing, there is a unique 
a c for which equality holds in fll5l) . Inserting the value (|) with q = 2 for a c 
yields — together with the second requirement on Zx, in the relativistic 
bound — the claimed stability criterion. ■ 

4 Instability with the Free Dirac Operator 

In this section we shall discuss the Brown- Ravenhall model (0). That is 
to say we consider the energy expression ([!]) with zero vector potential in 
the definition of the allowed electronic states (^), i.e., we take only \l/ e 

N 

7^tv,o = A where 7i + is the positive spectral subspace of the operator 
— ol ■ zV + m/3. We shall prove that there is no stability is in this model by 
showing that for any (sufficiently large) particle number, N, and any a > 
the energy is unbounded below. In other words, "stability of the first kind" 
is violated. It is nevertheless true, however, that for any choice of particle 
numbers and nuclear charges there is always a sufficiently small, nonzero a 
such that the energy is bounded below by zero. 

Since the positive spectral subspace 7i + for the free Dirac operator is 
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not invariant under gauge transformations we see that this Brown-Ravenhall 
model is not gauge invariant. (The previous, modified model discussed in 
Section || is not only stable, it is also gauge invariant.) More precisely, the 
energy spectrum depends not only on B but in fact on the full gauge potential 
A. The Brown-Ravenhall model is therefore physically meaningfully defined 
only if we make a fixed choice of gauge. The natural choice is the Coulomb 
gauge (radiation gauge), 

V • A = 0, 

since in quantum electro-dynamics this gauge implies that electrons interact 
via the usual Coulomb potentials and the coupling to the transverse field is 
minimal, i.e., derivatives are replaced by covariant derivatives. 

The interesting quantity is the lowest energy that the system can have. 

Definition 1 (Energy) 

E N>K :=inf£o[*,*]- 

where the infimum is taken over all divergence free A fields, pairwise distinct 
nuclear locations Ri, and normalized, antisymmetric states $ G H.n,o- 
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4.1 Stability with Small a and Small Particle Number 

Since this result is not a main point of this paper we shall be brief — even 
sketchy. If a single particle \& is in the positive spectral subspace of -D(O) 
then the action of -D(O) on \I/ is the same as multiplication of each component 
by (p 2 + m 2 ) 1 / 2 in Fourier space. For such functions we see that D(0)ty) 
exceeds |V|\1/), so we may as well replace -D(O) by |V| and also drop the 
condition that \1/ belong to the positive spectral subspace of -D(O). 

The next step is to use the lower bound on V c in (^) so that the energy 
is now bounded below by a sum of one-body operators, in a manner similar 
to that in Section |3| with a c as in Theorem [[]): 



So[*,*]>(*,Y l (l--)\Vv\*) + >fc jA+— / B\ (16) 



(Note again that the 'effective spin' q is 2, as can be seen be repeating the 
above argument.) Here j(r) is the current in the state \I/ and it is trivially 
bounded above pointwise by the density p(r) in the state ^ (defined in Ap- 
pendix |B|) . Therefore the integral involving A is bounded below by 



N 
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Now j B 2 > j \VA\ 2 and this is not less than i^ 3 " 2 ||A|^ by Sobolev's inequal- 
ity where K% = 4 1 / 3 (37r) _1 / 2 7r -1 / 6 (see |TJ|, p. 367). Similarly, the kinetic en- 
ergy (\I>, ^2v=i * s bounded below by 1.63g -1 / 3 J p 4//3 , which was proved 
by Daubechies || and which follows from (||). If we use these inequalities 
and then minimize the energy with respect to the unknown quantity ||^4||6> 
we easily find that the energy is non-negative as long as 



1.63(1 - a/a c ) > 2^N 2 / 3 Klq 1 ' 3 



a 



with q = 2. 

We shall show in Subsection |4.2| that the condition that N 2 ^ 3 a is small, 



which — as we just proved — ensures boundedness from below, is in fact also 
necessary for the energy to be bounded from below. 

4.2 Instability for All a and Large Particle Number 

The main result of this section is that there is no stability in this model for 
any fixed, positive a if N and K are allowed to be arbitrary. 

Theorem 2 (Instability) There exists a universal number C > such that 
for all values of the parameters a > 0, m > 0, K — 1, 2, 3, . . . ; and all values 
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of N = 1,2,3,.. and of Zi, Z 2 , . . . , Z K satisfying 



K 



K 



J2 Z *> Cmzx{a- 3/2 ,l 



}, iV>Cm a x{a- 3 / 2 ,l}, J2 Z > 2 



K=l 



we have 



E n ,k = -oo. 



Proof. The theorem follows if, for all E > 0, we show the existence of 
three quantities for which £ [&, < —E, with $ = ^ A • ■ ■ A -^at: 

A. A vector potential A with V • A = 0. 

B. Orthonormal spinors ipi, . . . , ip^ G 7i + . 

C. Nuclear coordinates Ri, . . . , R^. 

Our construction will depend on four parameters (to be specified at the 
very end), 5 > a momentum scale, which we shall let tend to infinity, 
6 > 0, which will be chosen sufficiently small (but independently of N), and 
P, Aq > which will be chosen as functions of N. Finally we denote by 
n 1: n 2 , n 3 the coordinate vectors (1, 0, 0), (0, 1, 0), (0, 0, 1) respectively. We 
shall use the notation that u> p = p/|p| is the unit vector in the direction 



p e M 3 . 



23 



A. The vector potential. We choose the vector potential A to have Fourier 



transform 



A(p) := A)Xb(o,55) (p) ( n 2 • w p ) n 3 x w p , 



where X-B(o,5<5) denotes the characteristic function (in p-space) of the ball 
B(0, 55) centered at with radius 55. Note first that A is real since A is 
real and A(p) = A(— p). Moreover, A is divergence free, i.e., it is in the 
Coulomb gauge, since — iV ■ A(p) = p ■ A(p) = 0. We easily estimate the 
self-energy of the magnetic field B = V x A corresponding to A 



-L |(V x A) 2 = -L J |p x A(p)| 2 rfp < Ap- 1 J%*d P = 2- 1 5 4 A 2 5 5 . 



Finally, we note for later use that 



(17) 



A(p) • nx = -A Xb(o,55)(p) (n 2 • w p ) 2 . 



B. The orthonormal spinors. For p G M 3 define 

/ 

Mpo (p) = ^/3/W)S- 3/2 



\ 



Xb(o,s)(p ~ Po) 




\ 
/ 



We then have a normalized ip Po G 7i + given by 

/ 

V^(p) = (2E(p)(E(p) + E(0)))- 1/2 



\ 



(E(p) + E(0))u po (p) 
p-au P0 (p) 



\ 
) 



(18) 



(19) 
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where E(j>) = (p 2 + m 2 ) 1 / 2 . Recall that this is the general form of a spinor 
in the positive spectral subspace 7i + for the free Dirac operator. 

For the sake of simplicity we shall henceforth assume that m — 0. We 
leave it to the interested reader to check the estimates for the general case 
m^O. We shall indeed consider spinors with momenta p such that we have 
p 2 (m 2 )~ 1 — > oo as 5 — > oo, i.e., E(p) « p. It is therefore straightforward 
to estimate the expressions in the general case m ^ by the corresponding 
expressions for m = 0. 

In particular, we have, for m = 0, 



^po(p) = 2- 1/2 



«Po(P) 



We shall choose iV points pi, . . . , G M 3 such that the following conditions 
are satisfied. 

1. min^ \p u - p M | > 25 

2. P < p v < 2P, for all v = 1, . . . , JV 

3. w p „ • ni > 1 - 9 2 , for all v = 1, . . . , iV 

Condition [I] ensures that the spinors ip Pl , . . . , ip PN are orthonormal. The 
importance of Conditions [| and will hopefully become clear below. 
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In order that Conditions and [I] are consistent with having N points 



(for large N) we must ensure that N balls of radius 5 can be packed into the 
domain defined by Condtions |] and Since small enough balls can fill at 
least half the volume of the given region we simply choose P such that 



Note that the assumption that N is larger than some universal number en- 
sures that the balls are small, i.e., that 5 is small enough compared to P. 
Thus, we have the condition 



In particular, since we shall choose 6 independently of N we may assume that 
N is large enough that the above condition implies P9 > 25. (Since we shall 
choose 5 — > oo we see that the momenta of the spinors satisfy p 2 m~ 2 — > oo.) 

We are now prepared to calculate J2u=i A/(A)\I/), where ^ = i/j Pi A 
■ ■ • A ip PN - We obtain 



2N < 



Vol ({p | P < p < 2P, 1 - 9 2 < lj p ■ m}) _ 7 2 P 3 
(4vr/3)(5 3 ~ 2 ~¥' 




(20) 




(21) 



where 



j v (r) := Vp„(r)aV P ,(r)) 



(22) 
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is the current of the u-th. one-electron state ip Pl/ , and 



T v :=ty Pv ,(-ia-V + /3m)1> Pv ) = j E(p)\u Pv (p)\ 2 dp < Vv + 5, (23) 



since we have assumed that m = and hence E(p) = p. 
We must evaluate the current integral 



jj u A = (2tt)- 3 / 2 2K 



KM ~ p)( a (p) ■ cr )( u; q • °")*v(q) rf p rf q 



(2tt)- 3 / 2 2K / / [A(q - p) • u; qWp >K„(q) 



+m£„(p) A(q-p) xcj q ) -(ru Pv (q) 



(24) 
dpdq. 



We first observe that 



J j iu * P M (A(q - P) x w q ) • <TW p „(q) rfprfq 



29fJ 



^p,(p) ( A (q - P) x u> q ) ■ n 2 a 2 M Pl/ (q 



dpdq = 0. 



The terms containing <j\ and vanish as they are clearly imaginary The 
term with o<i vanishes beacuse of the choice (|19D of u Pv . 

Note that u Vv (p)u Pv {c\) = unless |p — q| < 25 and |q — p u \ < 5. Thus 
w p „ - = P^iPu - q) + ^nPu 1 ^ - Vv) and we obtain for (p^ - q| < 25 
that 



|^ P ,-^ q |<2^; 1 <2(5P- 1 , 
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where we used that p v > P. Since ui Pi/ -n.i > 1—9 2 we have that |u> Pl/ — ni| < 9 
and hence 

A(q - p) • (« q - ni) < {25P- 1 + 9)A < 29A . 



Hence, since |p — q| < 25 we get from fll8|) that 



A(q - p) • u> q < -A [(u; q _ p ■ n 2 ) 2 - 29] 



Thus, 



Jj v ■ A = (2tt)- 3 / 2 2K //A(q -p) • w qUpi/ (pK„(q)dqdp 

[(u> q _ p • n 2 ) 2 - 26>] dpdq 



(2tt) 5 /2 



|q|.|p|«5 



^A 5 3 / / [( Wq _ p ■ n 2 ) 2 - 29] d P d q . 

W I J J | q |,|p|<l 



We now make the choice 



3 



and arrive at 



q|.|p|<i 



Idpdq 



lq|,|p|<l 



(w q -p ■ n 2 ) 2 dpdq 



Jj„-A<-^2n)-^A 
From ([H]) and (|23|) we therefore obtain 



N 



N 



(*,x;^(A)*)<x; 



i/=i 



|p,|+(5--(27r)- 1 / 2 Aoe ( 5 3 ^ 
o 



(25) 
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C. The nuclear coordinates. Finally, we show how to choose the nuclear 



coordinates following an idea in [18|. Consider the electronic density of the 



state V, p{r) = £f =1 |^(r)| 2 then 



A 



k=l " 1 1 «,A=1 



R 



A 



were we introduced D(p,p) := | f f p( r )| r — r'\~ 1 p(r')drdr'. 

Note now that J N~ x p = 1, i.e., N~ l p can be considered a probability 
distribution. We may therefore average (iff, V^) considered as a function of 
Ri, . . . , Rx with respect to the probability measure 

R 1; . . . , R K i ► iV-V(Ri) ■ ■ ■ N-V(Rjf). 

We obtain 

y K^)iV-V(Ri) • • • N-^R^dRi ■■■dR K 

K 

= [(l-(Z/N)r-N-^Z 2 K ]D(p,p). 



K=l 



We shall prove that [...] < 0. There are two cases. 

1) N > Z: By moving electrons to infinity we may assume that Z < N < 
Z + 1. Therefore 

K K 

[(l-(Z/N)) 2 -N- 2 Y,Zl] < Z- 2 -(Z + l)- 2 Y,Z 2 K < Z- 2 -2(Z + l)- 2 < 0. 

ft=l K=l 
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2) N < Z: We may move nuclei to infinity and assume that N < Z < 
N + max K Z K . Therefore 



K 



[(1 - {Z/N)) 2 - N~ 2 Z 2 K ] < (max Z K ) 2 N~ 2 - N~ 2 ^ Z 2 K < 0. 



K=l 



0. 



We can therefore find nuclear positions Ri, . . . , such that V c ^) < 



Using these coordinates together with ( |17D and (|25f) we get 



N 



p u + 5- -(2vr)- 1 / 2 A e^ 



•4iV\ 1/3 x . 



+ 2~ 1 5 4 A 2 <5 5 . 



We now choose P = (f^) 5 in accordance with ( p0|) . We then choose A 
such that |(27r)~ 1 /2^ oe5 3^ = 4 p ? Le>> 



At 



112V3 



— {2Ti) 1 ' 2 9-^e- l 5~ 2 N 1 '\ 



Then, using Condition ^| we have p u — (87i/3)Aoe5 3 6 < —2P. If we insert 
these values and estimates above we find 

£o[*, *] < 5 [-2N 4 / 3 (4/(70 2 )) 1/3 + iV + sVm- 2 / 3 ^- 10 / 3 ^ 1 ^/ 3 " 

It is now clear that the expression in [ ] will be negative for N sufficiently 
large. The energy can therefore be made arbitrarily negative by choosing 5 
large. | 
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A BKS Inequalities 

As a convenience to the reader we give a proof of some cases of the inequalities 
due to Birman, Koplienko, and Solomyak 0. The case needed in Section § 
corresponds to p = 2 below. There we are interested in (B — but here 
we treat (B — A) + to simplify keeping track of signs. The proof is the same. 
Recall that X + := (\X\ + X)/2. 

Theorem 3 Let p > 1 and suppose that A and B are two nonnegative, self- 
adjoint linear operators on a separable Hilbert space such that [B p — A p )^ p 
is trace class. Then [B — A) + is also trace class and 

tr(B - A)+ < tr(B p - 

Proof. Our proof will use essentially only two facts: X i— > X^ 1 is oper- 
ator monotone decreasing on the set of nonnegative self-adjoint operators 
(i.e., X > Y > =^ Y^ 1 > X^ 1 ) and X *—>■ X r is operator monotone 
increasing on the set of nonnegative self-adjoint operators for all < r < 1. 
Consequently, X i— > X~ r is operator monotone decreasing for < r < 1. 

As a preliminary remark, we can suppose that B > A. To see this, write 
B p = A p + D. If we replace B by [A p + then (B p - A p ) + = D + is 

unchanged, while X := B — A t— > [A p + D + ] 1 ^ p — A can only get bigger because 
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X i— > X 1//p is operator monotone on the set of positive operators. Since the 
trace is also operator monotone, we can therefore suppose that D = D + , i.e., 
B p = A p + C p with A,B,C> 0. Our goal is to prove that 

tr [{A p + C p ) 1/p - A] < tr C, (26) 

under the assumption that C is trace class. 

To prove (]26f) we consider the operator X := [A p + (7 p ] 1 / p — A, which is 
well defined on the domain of A. We assume, at first, that A p > e p for some 
positive number e. Then, by the functional calculus, and with 

E := [A p + C p ] {1 ~ p)/p and P := A x ~ p — E 

we have 

X = E[A P + C p ] - A l - p A p = -PA P + EC P (27) 

Clearly, P > and < P < e l ~P. 

Let Y := EC P . We claim that Y is trace class. This follows from 
Y*Y = C p E 2 C p < C p C 2 - 2p C p = C 2 . Thus, \Y\ < C, and hence tiY = 
ti C p l 2 EC p ' 2 < trC. 

It is also true that P is trace class. To see this, use the integral represen- 
tation, with suitable c > 0, A x ~ p = c f£°(t + A)~ 1 t 1 ~ p dt. Use this twice and 
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then use the resolvent formula. In this way we find that 

poo 

P = c (A p + t)- 1 C p (A p + C P + t)-H {1 - p)/p dt. 
Jo 

Since C is trace class, so is C p , and the integral converges because of our 
assumed lower bound on A. Thus, P is trace class and hence there is a com- 
plete, orthonormal family of vectors t>i, v%, . . ., each of which is an eigenvector 
of P. 

Since X > 0, the trace of X is well defined by Y^jLi( v ji Xv j) for any 
complete, orthonormal family. The same remark applies to EC P since it is 
trace class. Thus, to complete the proof of ( pop it suffices to prove that 
(vj,PA p Vj) > for each j. But this number is Xj(vj,A p Vj) > 0, where 
Xj is the (nonnegative) eigenvalue of P, and the positivity follows from the 
positivity of A. 

We now turn to the case of general A > 0. We can apply the above proof 
to the operator A + e for some positive number e. Thus we have 

tr [[(A + ef + C p ] 1/p -(A + e)] < tr C. (28) 

Let (fi, (p2, ... be an orthonormal basis chosen from the domain of A p . This 
basis then also belongs to the domain of A and the domain of \{A+e) p +C p ] l ^ p 
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for all e > 0. We then have 

tr X = ^(ifj.Xipj). 

3 

Note that a-priori we do not know that the trace is finite, but since the 
operator is non-negative this definition of the trace is meaningful. Operator 
monotonicity of X 1 ^ gives 

fo, [U + ef + Ct /P - {A + e)] tp,) > fo, (X - 

It therefore follows from (p8|) , followed by Fatou's Lemma applied to sums 
that 

trC > hminf^^^P + ^ + ^^-^ + e)]^) 
j 

> ^liminf(^, [[(A + e )* + - (A + e)] ip s ) 

3 

> tr X. 

■ 

B Counting Spin States 

Our goal here is to prove that when \1/ is in 7iN,Ai the antisymmetric tensor 
product of the positive energy subspace of the Dirac operator (with or with- 
out a magnetic field, A) then the one-body density matrix is bounded by 2 
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and not merely by 4, as would be the case if there were no restriction to the 
positive energy subspace. This result will allow us to use 2 instead of 4 in 
inequalities @ and (0). We thank Michael Loss for the idea of this proof. 

The one-body density matrix is defined in terms of an iV-body density 
matrix (or function) by the partial trace over N — 1 variables. We illustrate 
this for functions, but the proof works generally. If ^ is a function, then 



where z denotes a pair r, a and dz denotes integration over M. 3 and summation 
over the q 'spin' states of a. We are interested in q = 4, but that is immaterial 
for the definition. 

The kernel T is trace class; in fact its trace is qN. It is also obvi- 
ously positive definite as an operator. The first remark is that T < 1 
as an operator. To prove this easily, let ip be any normalized function 
of one space-spin variable z and define the function of N + 1 variables 
$(z Q , ...,z N ) := ^Oo^Oi, ...,z N ) + E^it-i^OjO^Oo, z N ), 
where Zj denotes the absence of Zj. This function $ is clearly antisymmetric 
and the integral over all variables of its square is surely nonnegative. How- 
ever, this integral is easily computed (using the normalization of ip and 




(N-l) 



^(r,cr, z 2 ,z 3 , . . .,z N )y(r', a',z 2 , • • • , z N )dz 2 ■ --dz N , 



to be (N + 1) - (N + 
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The next step is to consider the reduced kernel (without spin) defined by 

<i 

7(r,r') := ^ T(r, cr; r', cr), 
<j=i 

which evidently satisfies the operator inequality < 7 < q, since T < 1. 



The electron density referred to in Section 4.1 is defined by 



p(r) := 7(r,r), 

but it will not be needed in this Appendix. Another quantity of interest is 
the current, defined by 

j( r ) := 5^ r ( r ' (T;r ' r ) aCT ' r - 

<7,T 

It follows from this that |j(r)| < p(r) for every r G M. 3 . 

Our goal here is to prove the following fact about 7: 
// the N-body \1/ is in 1-Ln,a then the corresponding 7 satisfies < 7 < 2 as 
an operator. 

To prove this we introduce the unitary matrix in spin-space (related to 
the charge conjugation operator) 



/ \ 

1 



U 



\ 



-1 



/ 



where 1 denotes the unit 2x2 unit matrix. With a slight abuse of notation, we 
shall also use U to denote the U<S>1 acting on the full one-particle space, i.e., 
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(U)(r, a') = ^2 a U(a',a)f(r,a). The important point to note, and which 
is easily verified from the Dirac equation, is that ip G 7i + if and only if 
Uip G the negative spectral subspace of D(A). 

Given / G L 2 (R 3 ), we define F T to be the spinor F T (r,a) := f(r)5 a ^ T . 
Then evidently (/, 7/) = Y1 T (F T , TF T ). However, since the matrix U merely 
permutes the spin indices and possibly changes the sign from + to — , we 
have that £ T (F T , TF T ) = £ T (F T , IVF T ), with IV := U~ l TU. (Actually, 
the proof only requires that U be unitary, nothing more.) 

We claim that r + IV < 1 in which case we have proved that (/, 7/) < 
(7/2 = 2, as claimed. To see this, we note that T < 1 on 7i + and IV < 1 on 
Since the two subspaces are orthogonal, r + IV < 1 on the whole spinor 
space. 
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